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Perfect Space-Time Codes with Minimum and 
Non-Minimum Delay for Any Number of Antennas 

Petros Elia, B. A. Sethuraman and P. Vijay Kumar 


Abstract —Perfect space-time codes were first introduced by 
Oggier et. al. to be the space-time codes that have full rate, full 
diversity-gain, non-vanishing determinant for increasing spectral 
efficiency, uniform average transmitted energy per antenna and 
good shaping of the constellation. These defining conditions 
jointly correspond to optimality with respect to the Zheng-Tse 
D-MG tradeoff, independent of channel statistics, as well as to 
near optimality in maximizing mutual information. All the above 
traits endow the code with error performance that is currently 
unmatched. Yet perfect space-time codes have been constructed 
only for 2, 3,4 and 6 transmit antennas. We construct minimum 
and non-minimum delay perfect codes for all channel dimensions. 


I. Introduction 
A. Definition of Perfect Codes 

In [1, Definition 1], the concept of perfect codes is intro¬ 
duced to describe the n x n space-time codes that satisfy all 
of the following criteria: 

• Full rate. This corresponds to the ability of the code 
to transmit n symbols per channel use from a discrete 
constellation such as the QAM or the HEX constellation. 

• Full diversity. This corresponds to having all AXAX\ 
AX a difference matrix of the code, be non-singular. 

• Non vanishing determinant for increasing spectral effi¬ 
ciency. The determinant of any difference matrix, prior 
to SNR normalization, is lower bounded by a constant 
that is greater than zero and independent of the spectral 
efficiency. 

• Good shaping of the constellation. When the code is 

based on cyclic division algebras (CDA) ([8], [14], [1], 

[2]), the condition requires that the signalling set, in 

the form of the layer-by-layer vectorization of the code- 

2 2 

matrices, be isomorphic to QAM" or HEX" , where the 
isomorphism is given strictly by some unitary matrix. 

m Uniform average transmitted energy per antenna. The 
condition requires that the expected value of the transmit¬ 
ted power is the same for all antennas. In fact we will see 
that the structure of the code will allow for equal average 
power across the antennas as well as across time. 

The term perfect, coined in [1], was in reference to the 
ability of the codes to satisfy all the above criteria, as well as in 
reference to the codes having the best observed performance. 
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B. New results 

1) Additional properties satisfied by perfect codes: In ad¬ 
dition to the defining properties, perfect codes also satisfy 

• Approximate universality. This property was introduced 
in [18] to describe a code that is D-MG optimal [3] 
irrespective of channel statistics. Such codes exhibit high- 
SNR error performance that is given by the high-SNR 
approximation of the probability of outage of any given 
channel, and allow for the probability of decoding error 
given no outage to vanish faster than the probability 
of outage. Currently, the only family of approximately 
universal codes is that of cyclic-division algebras with 
non-vanishing determinant, the generalization of which 
was presented in [2]. Consequently perfect codes are also 
approximately universal. 

• D-MG optimality for any spacial correlation of the fading 
or the additive noise. This can be immediately concluded 
by post-multiplying the received signal matrix with the 
covariance matrix of the additive noise vectors and then 
using the approximate universality property of the code 
over channels with uncorrelated additive white noise. 

• Residual approximate universality. Any truncated code 
version resulting from deletion of rows maintains approx¬ 
imate universality over the corresponding channel. 

• Gaussian-like signalling. This is an empirical observation 
and it relates to a Gaussian-like signalling set with a 
covariance matrix that tends to maximize mutual infor¬ 
mation. 

• Information losslessness. This property relates to having 
full rate as well as unitary linear dispersion matrices [22], 
[20], and guarantees that the mutual information is not 
reduced as a result of the code’s structure. 

• Scalable sphere decoding complexity. This property guar¬ 
antees that as the number of receive antennas becomes 
smaller, the structure of the code allows for substantial 
reductions in sphere decoding complexity without essen¬ 
tial loss in performance. Eor MISO channels, the structure 
of the code will allow for reduction of sphere decoding 
complexity, from 0{'n?) to 0{n). 

2) Summary of presented contributions: In this paper we 
introduce explicit constructions of minimum-delay perfect 
space-time codes for any number n of transmit antennas 
and any number Ur of receive antennas. Non-minimum delay 
perfect codes are constructed for any delay T that is a multiple 
of n. We also generalize the defining conditions in [1] to be in¬ 
dependent of the channel topology, thus directly providing for 
joint maximization of the mutual information and approximate 


universality. For channels with a smaller number of receive 
antennas, this information theoretic approach allows for the 
construction of efficient variants of perfect codes which exhibit 
almost the same error performance as standard perfect codes 
but with substantially reduced sphere decoding complexity. 


II. Satisfying the Perfect-Code Conditions 

A. The general CDA structure 

Our constructions of the perfect n x n space-time block 
codes will be based on cyclic division algebras. 

As shown in other related works such as [ 8 ], [14], [1], 
[2], the basic elements of a CDA space-time code are the 
number fields F, L, with L a finite, cyclic Galois extension of 
F of degree n. For a being the generator of the Galois group 
Gal(L/F), we let z be some symbol that satisfies the relations 

£z = za{€) V £ G L and z^ = "f (1) 

for some ‘non-norm’ element 7 G F* := F\{0} such that the 
smallest positive integer t for which 7 * is the relative norm 
A’l/]P'(m) of some element u in L*, is n. The cyclic division 
algebra is then constructed as a right L space 

i:>=L©zL 0 ... 0 z"“^L. (2) 

A space-time code X can be associated to D by selecting 
the set of matrices corresponding to the representation by left 
multiplication on D of elements from a finite subset of D. 
For an arbitrary choice of scaled integral basis {/3i}"CQ^ of L 
over F (by a scaled integral basis we mean a set of numbers 
obtained by multiplying all the elements of an integral basis 
by the same nonzero real number), the n-tuple {fi,j}^~Q maps 
to 

n —1 

fbiPi. G L, G Ov (3) 

i=0 

where Op is the ring of integers of F. Consequently, prior 
to SNR normalization, the code-matrix X representing the 
division algebra element x = ^0 G L, is given 

by the defining equations (□3 to be 


£0 

7(t(£„_i) 

7cr^(£„_2) ■ 

•• 7f^”“H£i) 

£1 

o-(£o) 

7Cr^(£„_l) • 


£n-2 

Cr(£n-3) 


■■ 7a”-2(£„_i) 

£n-l 

Cr(£n-2) 

ct(£„-3) 

1 

0 


or equivalently 

n—1 

J =0 


diag{l ■ G) 


with 1. = [fj,o fj,i ■ ■ ■ where 


G 


/3o 


Pn-l 


a"-'(/3o) 


0 0 ••• 0 7 

1 0 ••• 00 

0 1 ••• 0 0 


0 0 ••• 10 


( 6 ) 


( 7 ) 


For this setup, it was shown in [14][9][2], that if 7 is inde¬ 
pendent of SNR and the fi j are from a discrete constellation 
such as QAM, then the code achieves 


1 ) full diversity 

2 ) full-rate 

3) non-vanishing determinant. 


1) Requirements for achieving equal power sharing and 
good constellation: Based on the above setup, we now show 
that using a unit-magnitude, algebraic, non-norm element 7 , 
and a unitary G, renders the code perfect. 

Let B := {/3o, • ■ • , Pn-i} be a basis for L/F such that the 


matrix 


G{B) = 


Po 


Pn-l 




is unitary. Furthermore, we assume that I 7 I = 1. 

Let F = Q(i) and let the {fij} be restricted to belong to 
the M^-QAM constellation: 


fij G -4qam = {a+ib I —(M—1) < a,b < M—1, a, 6 both odd}. 


For 


/ [/o ,0 fox ■ ■ ■ fo,n-l fl,0 ■ ■ • /l,rt-l ’ ’ ’ /n-l.n-l] 


we denote the code-matrix X in 0 as X{f) to emphasize 
that it is a function of the QAM vector /. If g is any function 
of the {fij}, we use E[p({/ij})] to denote the average 
77 ^ Then for ko G {0,1}, it is the case 

I-^QAmI f 

t«HQAM 

that 


E[|7'=V‘(£,)n 


i k 

i 

EE[|A,f|a‘(A)P] 

i 

i 

E[\hA\ 


( 5 ) It follows that the average energy per transmitted element of 
the code matrix is the same. 











With respect to the constellation shaping, let us now denote 
the layer-by-layer vectorization of X{f) as: 


vec(X(/)) := 


^nxn 

Onxn 


Onxn 


Onxn 

Onxn 


4 


0 


n X n 


.pl’T’-l) 


fn-1 


where 
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0 . 

0 
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1 

0 . 

. 0 

0 

pA) = 

0 

0 . 

1 

0 

r(2) = 

0 

0 . 

• 7 
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0 

0 . 

0 
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0 

0 . 

. 0 
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1 

0 

. 0 

0 
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. 0 
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ra-l) _ 

0 
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■ 7 
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We observe that each vector resulting from the layer-by-layer 
vectorization of any code-matrix, prior to SNR normalization, 
is exactly the linear transformation of the n^-tuple / from 
OAM”"" or HEX"\ by the unitary matrix 



G 

Onxn 

Onxn 

’ Onxn 


Onxn 


Onxn 

■ Onxn 

Rv = 

Onxn 

Onxn 

G • . 

' Onxn 


Onxn 

Onxn 

Onxn 



As a result the signalling set prior to SNR normalization is 
from the lattice 


Asignal = {/i?., /eQAM"'} (9) 

or Asignai = {/i?n, / G HEX”""}. Modifying / by hor¬ 
izontally stacking its real and imaginary parts, and 

respectively, and transforming the resulting vector / = 
A. LJ e by 



Rv,Re 

~Rv,Im 


Rv^Im 

Rv,Re 


( 10 ) 


we get the real and imaginary stacking of the signalling set to 
be from 


A_ai = {/ K, f e z"""}, = I, 


2n^X2n^ 


( 11 ) 


exactly satisfying the good constellation shaping condition as 
described in [1, Definition 1]. 


Equivalently, for any set of ki G {0,1}, 



2 

4 

2 

fo,j 

7 "^ ^(4) 

— 

o-(4) 

= 

h,J 







allowing for 

Tr(XtX) = ||vec(X(/))f = ||/f 
and making the collection 

{vec(X(/)) I / G 

represent a cubic constellation in -dimensional complex 
space that is isometric to 

We proceed to find a proper unit-magnitude non-norm 
element 7 , and then a proper unitary matrix G. 

B. Uniform average transmitted energy per antenna 

We now provide algebraic, unit-magnitude ‘non-norm’ ele¬ 
ments 7 for suitable cyclic Galois extensions L/F, independent 
of SNR. We will henceforth denote the primitive root of 
unity as toi, i.e. w/ = and by k* we will denote the 

complex conjugate of fc G C. We directly state the construction 
method for the different cases of interest, depending on 
whether the base field is F = Q(j) (QAM) or F = Q(a; 3 ) 
(HEX). 

Proposition 1: {Construction of the non-norm element for 
the QAM code) Let n = 2®ni where ni is odd. Then there 
exists a prime p congruent to 1 mod rii. Furthermore, there 
exists a prime q that is congruent to 1 mod 4, as well as 
congruent to 5 mod 2®+^, and which has order ord(g) = ni 
and splits in Z[z] as g = ttitt} for a suitable prime tti G Z[z]. 
The fields Q(a;p) and Q(i) are linearly disjoint over Q. Let 
K be the unique subfield of Q{i){ujp) of degree ni over Q(z) 
and let L = K • Q(a; 2 b+ 2 ). Then L is a cyclic extension of 
Q(z), and the element 


TTl 



is an (algebraic) unit-magnitude element that is a non-norm 
element for the extension L/Q(i) and is independent of SNR. 

1 2i 

(When ni = 1, we take L = Q(a; 2 b+ 2 ) and 7 =-.) The 

n X n matrices arising from equations i} and ^ with 
the above choice of F, L and 7 , and with any choice of scaled 
integral basis {f3i}^ZQ hence yield a full-diversity, full-rate 
code over QAM with non-vanishing determinant satisfying the 
additional equal power sharing constraint. 

Proof: See Appendix |I] ■ 

Proposition 2: {Construction of the non-norm element for 
the HEX code) Let n = 2®ni, s G {0,1}, where ni is odd. 
Then there exists a prime p > 3 congruent to 1 mod ni. 
Furthermore, there exists a prime q that is congruent to 1 mod 
3 and which has order ord(g)| 2 * = ni and splits in Z[a; 3 ] 
as g = ttitt} for a suitable prime tti G Z[w 3]. If s = 1 then 
g should also be congruent to 3 mod 4. The fields Q(a>p) 
and Q(a; 3 ) are linearly disjoint over Q. Let K be the unique 





























sub field of Q(a; 3 )(a;p) of degree ni over Q(aJ 3 ) and let L = 
K • Q(a' 2 <!+i). Then L is a cyclic extension of Q(a; 3 ), and the 
element 

TTl 


7 = 


is an (algebraic) unit-magnitude element that is a non-norm 
element for the extension L/Q(a; 3 ) and is independent of 
SNR. (When ni = 1, so s = 1, we take L = Q(a; 3 )(^), 

and 7 = —i—|.) The n x n matrices arising from equations 
3 

( 0 , 0 and 0 with the above choice of F, L and 7 , and 
with any choice of scaled integral basis {/3i}"7g^ hence yield 
a full-diversity, full-rate code over HEX with non-vanishing 
determinant satisfying the additional equal power sharing 
constraint. 

Proof: See Appendix U ■ 

Consequently, we have constructed algebraic, unit- 
magnitude ‘non-norm’ elements, valid for use in perfect 
codes, for any n. Some examples are given in Table |I] 


TABLE I 

Non-Norm Elements 


No. of Antennas 

Non-norm ‘ 7 ’ 

2 

{2 + i)/(l+2i) 

3 

(3 + i^s)/ ( 3 + CQg ) 

(1 + 9g; 3)/(9 + cJs) 

4 

(2 + 1)/{2-1) 

5 

(3-|-2i)/(3 - 2i) 

6 

(3 + 7g;3)/ (3 + 7uj^) 

7 

(8 -H 5i)/{5 + 8i) 

8 

(2T7)7{rT2l) 

9 

(3 -1- ivsi/ll + 3 uJ3) 

(4 9i)/(9 -H ii) 


C. Good Constellation shaping 

For any dimension n, we now proceed to describe the con¬ 
struction of the unitary matrix G that complies with the cyclic 
Galois requirements of the division algebra. The construction 
method will involve the embedding of the scaled integral basis 
of a submodule of Ol over the ring of integers Of of F, where 
F,L are as in Propositions ^ and |2] 

We proceed to first construct lattices for any odd dimension 
ni, then lattices of dimension 2 ®, s G and then proceed 
to combine them in order to give the final desired lattices for 
any dimension n over Q(z) or any dimension that is not a 
multiple of 4 over Q(u; 3 ). Without any loss of generality, we 
will be analyzing the QAM case, corresponding to F = Q(^). 
Unless we state otherwise, the same results will also hold for 
F = 

I) Orthogonal Lattices in a Cyclic Galois Extension over 
Q of Odd Degree : Recently, the authors in [12], Section 
V, give a detailed exposition of a previous result in [10] of 
an explicit construction of 5 -dimensional orthogonal lattices 
that belong in a 5 -degree cyclic Galois extension K' over Q, 
with the restriction that 5 be an odd prime integer. We here 
show that the same construction actually gives ni-dimensional 
orthogonal lattices in O'^, for any odd integer ni. Moreover, 
the held K'(i) will be precisely the held K of Proposition Q] 


The steps for constructing an ni-dimensional orthogonal 
lattice in Ok are as follows: 

« pick a guaranteed to exist odd prime p = 1 (mod ni) 

277% 

• let w = Wp = e p 

• hnd a guaranteed to exist primitive element r of Z* 

« for m = create a = ~ where 

rP-i = 1 

• Find a guaranteed to exist A such that A(r — 1) = 1 (mod 
p) and let 2 : = uj^a{l — uj) 

P-1 

. For cr(w) = wU let x = J2k=i (^)- 
The element x is hence in the held K', the subheld of Q(uj) 
hxed by (t"i, of degree ni over Q. It is then the case that 
the following lattice generator matrix Gm is unitary and the 
resulting lattice (which arises from the canonical embedding 
of the free Z module generated by cr(x)/y/p, ..., 

in K”) is orthogonal. 



X 

cr(x) • 

■■ a"i- 2 (x) 


G 

^ni — 

p 

cr(x) 

cr^x) • 

•• cr”i-l(x) 

X 

cr^x) 

aHx) ■ 

X 

cr(x) 


cr”i-l(x) 

X 

•• cr”i-7x) 

cr”i-7 


Since Q{oj) and Q(z) are linearly disjoint over Q, the held 
K = K'(z) will be cyclic over Q(t), and the elements xj^/p, 
a{x)jy/p, ..., (T"i“^(x)/y/p will be a scaled integral basis 
for K/Q{t). 

Proof: See Appendix HIl ■ 

More specihcally the hrst row of Gm is given by 

p-1 

P 

and the rest of the circulant matrix by: 

G„i(i-I-l,j) = + 1 mod ni), i = 0, • • • , ni - 2. 


Example 1: The hrst row of the 9-dimensional Gg is 
(-2.831 7.298 - 1.435 4.149 - 8.688 - 8.451 - 

6.414 5.355 — 7.983) and every next row is obtained by 
a single left cyclic shift of the previous row. The matrix 
was obtained by setting rii = 9, p = 19, r = 3 and 
A = 10. Similarly the hrst row of the 15-dimensional G 15 
is ^(-2.242 6.361 - 10.78 -8.071 7.253 -9.45 1.127 - 

3.334 8.806 -4.391 10.442 5.404 -11.12 -11.004 -9.989 ) 
obtained by setting ni = 15, p = 31, r = 3 and A = 16. 

2) Lattices of dimension m = 2® [13]: For when the 
information set is QAM, then F = Q(i) and we consider 
K. = Q(ujm) where M = 2®+^ and lum = w = the 

primitive root of unity. Q(a;) is a cyclic Galois extension 
over Q(i). Considering that the order of 5 in Z]^=GaZ(K/Q) 
is m = 2® = we see that for a G Gal{Q{io)/Q) 

such that a{Lu) = w®, it is the case that a{i) = a{uf ) = 

= I which gives that 
Ga((K/Q(i)) =< cr >. Taking to 

















be the integral basis over Q(j), the canonical embedding then 
gives the lattice generator matrix 


Ge 




1 






(13) 


The fact that the lattice corresponds to the ring of integers 
of the m-dimensional cyclic Galois extension K over Q(t), 
allows for Ge to be directly used in to construct the mxm 
space-time code. 

Now for Ti = [1 u! 

0 , 1 , • • • ,m — 1 , being the i 


,5‘2 


,5’3 


have that — 

M 




... ijj- 

row of 


= Since 

5 has order ^ in then 5* ^ 5-^ Vi ^ j, i, j = 

0,1, • • • , ^ — 1. This combines with the fact that fc(5* — 5-^) = 
k5^{5^~^ — 1) = 0 (mod 4) so that each summand pairs with 
another summand in the summation so that their ratio is 
This symmetry, the fact that -^ = 0 (mod 4) and the fact that 
uj^ + ) "S" =0, means that each summand oj^ has another 

summand as its additive inverse. Together with the fact that 
the complex conjugate of oj is results in = mSij 
and in the desired orthogonality GeG\ =■ I. The lattices apply 
only for codes over QAM. 

3) Combining lattices: We will need the following, which 
is an easy modihcation of Proposition 6 in [12] and which 
eventually guarantees for the creation of lattices over a cyclic 
Galois extension for any dimension n over Q(^), and any 
dimension that is not a multiple of 4 over <Q{uj 3 ). 

Lemma 3: Let L be the compositum of I Galois extensions 
Ki over Q of co-prime degrees n^. Assuming that there 
exists an orthogonal Ok .-lattice generator matrix Gi for all 
i = 1, 2, • • • , I then the Kronecker product of these matrices 
is a unitary generator matrix of an n-dimensional lattice in 
Ol, n = 

For this, the discriminants are not required to be coprime 
since the involved helds already have coprime degrees, so 
their composite is their tensor product over Q. Specihcally, for 
F = Q(i), for any n = ni2®, ni odd, the orthogonal lattice 
generator matrix G is the Kronecker product of the generator 
matrix of the ni-dimensional lattice from Section Hl-C. 1 1 and 
that of the cyclotomic lattice of dimension 2®. For F = Q(w 3 ), 
for n = ni odd we again use the rii-dimensional lattice from 
Section iri-C.ll and for n = 2ni, ni odd, the orthogonal lattice 
generator matrix G is the Kronecker product of the generator 
matrix of the ni-dimensional lattice from Section in -C. II and 

'1 i ' 

1 —i 

For F = Q(a; 3 ), for n = ni odd we again use the ni- 
dimensional lattice from Section Hl-C. II and for n = 2ni, ni 
odd, the orthogonal lattice generator matrix G is the Kronecker 
product of the generator matrix of the m-dimensional lattice 

from Section liTO] and the matrix G 2 = 


i = 

in disk we 


jk—0 

4>{M) 

2 


matrix G 2 = 


1 i 

1 —i 


coming 


from the held Q(?). 

The above orthogonal lattice generator matrices correspond 
to a suitable scaled integral basis of the n-dimensional cyclic 
Galois extension L/F, dehned (respectively) in Propositions 
n and I 2 I As discussed above, these matrices allow for good 
constellation shaping. Consequently, with this choice of lattice 


generator matrix and the choice of L/F and 7 as in Propo¬ 
sitions [H and I 2 I the code dehned by equations (l3i7> form 
perfect codes satisfying full-diversity, full-rate, non-vanishing 
determinant, equal power sharing, and good constellation 
shaping. 

111. Information theoretic interpretation and 

GENERALIZATION OE THE PERFECT CODE CONDITIONS 

The D-MG tradeoff [3] bounds the optimal performance of 
a space-time code X operating at rate R bpcu, corresponding 
to a multiplexing gain 

R 

log 2 (SNR)- 

The diversity gain corresponding to a given r, is dehned by 

Hr) = - lim . 

^ ^ SNR^oo log(SNR) ’ 

where Pe denotes the probability of codeword error. For the 
Rayleigh fading channel, Zheng and Tse [3] described the 
optimal tradeoff between these two gains by showing that for 
a hxed integer multiplexing gain r, the maximum achievable 
diversity gain is 

d[r) = (n — r){nr — r). (14) 

The function for non-integral values is obtained through 
straight-line interpolation. 

We use this D-MG approach as a basis for interpreting and 
generalizing the conditions that dehne perfect-codes. 

4) Full rate condition: Consider an n x T code X where 
each code-matrix carries m information symbols per channel 
use from a discrete constellation A such as QAM. It is then 
the case that 

lA”! = 2^^ = 2 ’'^'° 82 Snr ^ sNR’'^ = \A\^'^ 

which implies that |Vl|=SNR™ and since the constellation 
is discrete, we have that E[||a G Vl||^] = |Vl|. The fact that 
each element X^ j of a code matrix is a linear combination of 
elements of A, gives that 

E[||X,,,f] = |^|=SNR^. 

The SNR normalizing factor 1 / that guarantees that 
E[||i/iFX|j|.] = E[||i/Ar|||.]=SNR is then given by 

i/2=SNR1-™. (15) 

Without loss of generality we can assume that there exist 
two code-matrices Xi,X 2 G X, with each Xi mapping the 
information nm-tuple {ai,0, O,--- ,0}, where at = SNR°. 
As a result, the determinant and trace of the difference matrix 
AX, is a polynomial of degree less than n over a = ai — 
a 2 = SNR°, with coefficients independent of SNR, i.e. 

det(AXAX^) = Tr(AXAXl') = SNR° 

and thus with all its eigenvalues 

h = SNR°. 























The corresponding pairwise error probability PEP(Xi ^ X 2 ), 
in the Rayleigh fading channel [4], [7], then serves as a lower 
bound to the codeword error probability Pg, i.e., 


Pe > PEP(Xi^X2) 


n 





which results in a diversity gain of 

d(r) < rirU (1 -I . 

V TO/ 

What this means is that to discrete information symbols per 
channel use can potentially sustain reliable communication for 
up to rate i?max ~ TOlog 2 (SNR). 

Eor large SNR, the outage capacity over an n x Rayleigh 
fading channel is given by Cout ~ min{n, rir} log 2 (SNR), 
implying a maximum achievable multiplexing gain of 
rmax = min{n, rir}. Consequently the relation between i?max 
and Cout, allows for the interpretation that the full rate defin¬ 
ing condition is necessary for reliable transmission at rates 
close to the outage capacity of the Rayleigh fading channel, 
independent of the channel topology. Equivalently, given some 
rate R, the full rate defining condition is necessary for reliable 
transmission at the smallest allowable SNR 

SNRmin = 

again independent of the channel topology. Let us now re¬ 
examine the full-rate condition, in conjunction with the deter¬ 
minant condition. 

5) Non-vanishing determinant condition: We consider the 
n X T truncated code X, T > n, constructed by deleting 
the same T — n rows from all the code-matrices of a 
T X T perfect code X . We have seen that for an n x n code 
mapping information symbols from a discrete constellation 
{n information symbols per channel use), the standard n- 
dimensional ‘folding’ {\X\ = |-4|" ) forces a normalizing 
factor of = SNR^“", whereas in the truncated nxT code 
mapping information symbols information symbols 

per channel use), the constellation is folded in T dimensions 
(\X\ = requiring for 

z/2 = SNR^"T. 

This scenario accentuates the fact that in essence, we are 
limited by a lower bound on the determinant of the energy- 
normalized difference matrix AXAXK As a result, for the 
n-dimensional case, the defining condition of non-vanishing 
determinant for the non-normalized matrix > SNR°, 

translates to 

det[z/^AAAAt] > (i/2)”SNR° = (SNR^"-)" 

= SNR'""’' 

which, for the nxT case with T-dimensional folding, trans¬ 
lates back to the determinant bound 

det(AAAAl') > SNR-’r('^-”) 

for the non-normalized code-matrices. But from [2], [16], 
we see that the above determinant bound is the best that 


any code can attain, thus allowing us to generalize the full- 
rate, full-diversity and non-vanishing determinant perfect code 
conditions, to the general condition of having 

de%[v'^AXAX'^\ > SNR'""’’, 0 < r < min(n,ni.). (16) 

In regards to non-minimum delay perfect codes, let us 
briefly note that codes resulting from row deletion of perfect 
codes essentially maintain all the conditions of the original 
minimum-delay perfect code constructions except that now the 
vectorization of the code-matrices is not isometric to QAM^ . 
Non-minimum delay perfect codes can be constructed though 
for delays T = nk, k G that are multiples of n, by 
the horizontal stacking construction found in [2], [16] which 



the isometry of the code matrices with QAM" . 

Let us now incorporate all the perfect-code defining condi¬ 
tions in order to provide an information theoretic interpretation 
that spans both the high and the low SNR regimes. 

6) Approximate universality, information losslessness and 
Gaussian-like signalling: We begin with: 

Theorem 4: Perfect codes are both approximately universal 
as well as information lossless. 

Proof: See Appendix EYl ■ 

The code’s information losslessness, shown in the proof to be 
the result of the CDA structure and the last two conditions, es¬ 
sentially allows for the code to maintain the maximum mutual 
information corresponding to the channel and signalling set 
statistics. This mutual information is empirically related to the 
Gaussian-like signalling set and its good covariance properties, 
observed in Eigure ^ The expedited rate with which the 




Fig. 1. Gaussian nature of the signalling set of the 3x3 perfect code, 
compared to QAM and random Gaussian signalling (left). Covariance of 
columns of perfect codes (center) compared to covariance of random Gaussian 
vectors (right). 


signalling becomes Gaussian, relates to the high-dimensional 
and orthogonal nature of the lattice generator matrix which 
together with a unit-magnitude non-norm element, jointly 
allow equal magnitudes for the diagonal elements of the 
covariance matrix of the signalling set. 

Let us now draw from the information theoretic interpreta¬ 
tion of the defining conditions and provide variants of perfect 
codes that are specifically tailored for channels with a smaller 
number of receive antennas, and which manage to maintain 
good performance at a considerably reduced sphere decoding 
complexity. 

A. Channel topology and efficient variants of perfect codes 

We have seen that nxn perfect codes utilize n different lay¬ 
ers to achieve approximate universality for all n^. Each layer 
has non-vanishing product distance and maps n-elements from 













a discrete constellation, thus maintaining two properties that 
were shown in [18, Theorem 4.1] to guarantee for optimality 
over the statistically symmetric parallel channel, i.e. a channel 
with a diagonal fading coefficient matrix, as well as potentially 
allowing for optimality over the statistically symmetric n x 1 
MISO channel. To offer intuition, we observe that the sum- 
capacity of the rir independent MISO channels relates to the 
full rate condition, whereas the achieved full diversity relates 
to the CDA structure and the discreteness of the powers of T 
which manage to translate the non-vanishing product distance 
to an overall non-vanishing determinant, and thus to keep 
the different layers independent and at some non-vanishing 
distance from each other. The full rate condition comes with 
a sphere decoding complexity of O(n^), but as the number of 
MISO channels reduces with Ur, so does the required decoding 
complexity. Codes over such channels can have the form 

n—1 

X = 

j=Q 

where T{i,j) = I, i = j e [l,..,n^] and T{i,j) = 
0 otherwise, or can have the form 

j=o ^ ^ 

Note here that the above codes have not been proven to be 
D-MG optimal. 

Motivated by the down-link requirements and by the cooper¬ 
ative diversity uses of space-time coding in wireless networks, 
we will concentrate on the MISO case {n^ = 1), for which a 
D-MG optimal perfect code variant 

Xd = {diag{^ = diag{[_ ■ G), V/ e QAM"}. (17) 

with sphere decoding complexity of 0{n) was recently con¬ 
structed in [49], [50] for all n, together with the code 

n — 1 

X„ = {X = Y, fkT\ V/fc e QAM-HEX} (18) 

k=0 

that corresponds to the center of the division algebra. With 
the exception of n = 2, Xi^. has not yet been proven to be 
D-MG optimal. Figure |2] provides a performance comparison 
between the single dimensional perfect code variant with the 
equivalent standard perfect-code. 

IV. Recent developments involving pereect codes 
The proposed high-dimensional perfect codes have had an 
impact on establishing outage-based optimality expressions 
for wireless networks where independently distributed users 
utilize different parts of space-time schemes to relay messages 
for one another, hence improving the overall quality of service 
([43] etc). Up to now, outage-based optimality results were 
known only for infinite time duration networks in which 
the assisting relays required full knowledge of the channel. 
Encoding was based on random Gaussian codes. Using perfect 
codes as an information theoretic tool, it was shown in [47], 
[49] that the same optimality can be achieved, for finite 
and minimum delay, and without requiring knowledge of the 


diag{l -T-G) 



Fig. 2. For n = 2,nr = 1, the single-layer perfect code variant exhibits 
similar error performance as the 2-layer (standard) perfect code. This changes 
when a second receive antenna is added. 


channel at the intermediate relays. This was achieved for the 
most general network topology and statistical characterization. 

Perfect codes and some perfect-code variants, provided for 
the first ever optimal encoding method [47], [48], [49], [50] 
in several cooperative-diversity schemes such as the non¬ 
dynamic linear-processing (receive-and-forward) scheme [44], 
the non-dynamic selection-decode-and-forward scheme [45] 
and finally for the dynamic receive-and-forward scheme [46]. 


V. Examples of new perfect codes and simulations 


A. Examples of new perfect codes 

• A 2 X 2 perfect code can be chosen to have code-matrices 
which prior to SNR normalization, are of the form 


X 



/o ,0 + /OA^S 

/l,0 + /l.lWg 

/o,o + /o,iw| 
/l,0 + /l,iw| 


7(/i.o + /i,icr(u;i)) 
/o ,0 + fo,lO-{ujg) 

7(/i.o + /i.it4) 
fo,0 + fo,1^8 


where fij are from the desired QAM constellation, ws := 
and 7 = Matrices map nf — 4 information elements 

from QAM. Furthermore the signalling set, in the form of the 
layer-by-layer vectorization of the code-matrices, before SNR 
normalization, comes from the lattice 


^ = {[/o,0 fo,l fl,0 fl,l]Rv ■ 

''^[/o.O)/o.ij/i.o,/i,i] G QAM" } 

where 

110 0 
a;| Wg 0 0 

0 0 1 7 
0 0 Wg 7a;g 

satisfying the defining condition of 
RvRI = h- 

We find the smallest possible determinant, prior to SNR 
normalization, to be 

det(AXAXt)^i„ = 1 














which is larger than some previously constructed 2 x 2 
perfect codes. The code’s performance improves if the existing 


G = 


1 1 
,3 


is substituted with G 2 = 


0.5257 0.8507 

0.8507 -0.5257 


, ^8 

taken from [12 
Other examples: 

. The 5 X 5 perfect space-time code is given by 


X=\x = J2^^{dtag{l^-G,)), 

3=0 


£. G QAM^ 


for r given in 0 based on 7 = , and generator matrix 


- 0.3260 

0.5485 

- 0.4557 

- 0.5969 

- 0.1699 

0.5485 

- 0.4557 

- 0.5969 

- 0.1699 

- 0.3260 

- 0.4557 

- 0.5969 

- 0.1699 

- 0.3260 

0.5485 

- 0.5969 

- 0.1699 

- 0.3260 

0.5485 

- 0.4557 

- 0.1699 

- 0.3260 

0.5485 

- 0.4557 

- 0.5969 


. The 7 X 7 perfect space-time code is given by 

r ® 

X = lx = {diagil^ • Gr)), 4. G QAM" 

^ j=o 

for r based on 7 = and generator matrix G 7 = 


- 0.681 

0.163 

- 0.449 

0.077 

0.082 

0.276 

- 0.469 

0.163 

- 0.449 

0.077 

0.082 

0.276 

- 0.469 

- 0.681 

- 0.449 

0.077 

0.082 

0.276 

- 0.469 

- 0.681 

0.163 

0.077 

0.082 

0.276 

- 0.469 

- 0.681 

0.163 

- 0.449 

0.082 

0.276 

- 0.469 

- 0.681 

0.163 

- 0.449 

0.077 

0.276 

- 0.469 

- 0.681 

0.163 

- 0.449 

0.077 

0.082 

- 0.469 

- 0.681 

0.163 

- 0.449 

0.077 

0.082 

0.276 


> The 25 X 25 integral restriction code 
is given by 


r 

Xir = lx = SfeGQAM 

fc =0 

with 7 = This code has the same sphere decoding 

complexity of 0(25) as the 5x5 standard perfect code in the 
example above, and is expected to have the same performance, 
when Ur = 1, as the 25 x 25 perfect code whose sphere 
decoding complexity is 0(625). 


B. Simulations 

All the simulations assume CA/’(0,1) fading and thermal 
noise. A sphere decoder was used. We begin with Figure |3 to 
indicate the performance improvement as the different dehning 
conditions are satisfied one-by-one. The first curve from 
the top corresponds to satisfying the full-diversity condition 
(commutative CDA code - orthogonal design). The second 
curve now includes the full-rate condition (random, full-rate, 
linear-dispersion codes). The third curve corresponds to the 
family of D-MG optimal but not information lossless CDA 
codes presented in [ 2 ], which achieve the hrst three criteria 
of full-diversity, full-rate, and non-vanishing determinant. The 
performance transition from the CDA codes to perfect codes 
is described by the next two curves. Figures 0 and |3 show a 
comparison of the 2 x 2 unified perfect code presented here, 
with some perfect codes from [1] and with the Alamouti code 
(rir = 2). The Golden code [15] performs best among all 
existing 2x2 perfect codes. When rates are lower, the unified 
perfect and the Golden code perform better than the orthogonal 



Fig. 3. Performance improvements attributed to achieving the different 
criteria for the perfect codes 


design whereas one of the perfect codes does not always 
do so. For higher rates, all considered perfect codes perform 
substantially better than the orthogonal design. At all rates and 
all SNR, the perfect code constructed here has performance 
very close to that of the Golden code. In Figure]^ we show the 
performance of the newly constructed 5-dimensional perfect 
code and compare that with the corresponding 5x5 single¬ 
dimensional commutative perfect code (HI. As expected, the 
former utilizes fully the = n = 5 channel and is thus able 
to transmit with a small probability of error at high rates and 
low SNR. 



Fig. 4. Lo\v rate compaiison of the unified perfect code with two perfect 
codes from [1] and with the Alamouti code 


VI. Conclusion 

We have explicitly constructed perfect space-time codes for 
any number n of transmit antennas, any number Ur of receive 
antennas and any delay T that is a multiple of n. Achieving all 
the defining conditions from [ 1 ], allows for perfect codes to 
exhibit performance that is currently unmatched. The informa¬ 
tion theoretic interpretation of the exhibited good performance 
both for low and high SNR, is that the defining conditions 
jointly endow the code with approximate universality and the 
ability to provide for near optimal mutual information. 

High dimensional perfect codes cover a much needed re¬ 
quirement for optimal codes in multi-user cooperative diversity 
wireless networks, where each user acts as a transmit antenna. 
Specifically, perfect codes have already been used to establish 
























Fig. 5. High rate comparison of the unified perfect code with two perfect 
codes from [1] and with the Alamouti code 



Fig. 6. Compaiison of the 5x5 perfect code with the 5x5 single dimensional 
integral restriction code (nr = 5). 


the high-SNR outage region of unknown channels, and have 
provided the first ever optimal schemes for a plethora of 
cooperative diversity methods. 

Appendix I 

Proof of construction methodology for 

NON-NORM ELEMENTS 

We will prove here Propositions 11121 

For future reference, we first recall three results that relate 
to identifying a “non-norm” element 7, i.e, an element 7 G F* 
satisfying 7® ^ 0 < i < n for some n-dimensional 

field extension L of F. 

Lemma 5: [9] Let L be a degree n Galois extension of a 
number field F and let p be a prime ideal in the ring Of below 
the prime ideal ^ C Of with norm given by ||^p|| = Ijplj-^, 
where / is the inertial degree of ^ over p. If 7 is any element 
of p \ p^ then 7® ^ iVL/F(L) for any i = 1, 2, • • • , / - 1. 

In order to find a “non-norm” element 7 in F = Q(z) (F = 
Q(w 3 )), it is sufficient to find a prime ideal in Z[i] (Z[w 3 ]) 
whose inertial degree / in L/F is / = [L : F] = n. Such an 
ideal is said to be inert in L/F. 

Lemma 6: [38] Let p be any odd prime. Then for any k G 
Z, Z* is cyclic of order 4){p^). For any integer / dividing 
there exists an a G Z*^ such that a has order / in Z*^. 

Theorem 7: (Dirichlet’s theorem) Let a, m be integers such 
that 1 < a < m, gcd(a, m) = 1. Then the progression 


{a, a -I- TO, a -I- 2to, ..., a -I- km,.. .} contains infinitely many 
primes. 

We will now proceed to establish the exact methodology 
that will give unit-magnitude non-norm elements 7 , for the 
different cases of interest. 

a) Unit-magnitude, non-norm elements for F = Q(i): 
Let 

r 

n = 2 ® n pf = 2 ®ni (say) 

i=l 

where pi are distinct odd primes. Assume first that ni > 1. 
Let p be the smallest odd prime p such that ni \ (p — 1). 
The cyclic group Z* contains an element whose order equals 
(p — 1). Let a denote this element. Our first goal is to find a 
prime q such that 

q = 5 (mod 2 *+^) 

q = a (mod p). 

Note that 

<7 = 1 (mod 4) . 

Since (2®+^,p) = 1, we can, by the Chinese Remainder 
Theorem, find an integer b such that 

b = 5 (mod 2"+2) 
b = a (modp). 

Note that such an integer b is relatively prime to 2^'^^p. 
Consider the arithmetic progression 

b + l{2^+^p), I = 0,1,2,... 

By Dirichlet’s theorem, this arithmetic progression is guaran¬ 
teed to contain a prime q having the desired properties. Now 
let us verify that this leads to a CDA. 

Let K be the subfield of Q(u;p) that is a cyclic extension 
of Q of degree ni. Let K be the compositum of IK and Q(*) 
and let L be the compositum of the fields K and Q(a; 2 <i+ 2 ). 
Note that L is cyclic over Q(i), since it is a composite of the 
cyclic extension Q(u; 2 » + 2 )/Q(i) of degree 2 ®* and the cyclic 
extension K/Q(i) of ni (note that 2 ® and ni are relatively 
prime). Now consider the decomposition of the prime ideal 
(q) in the extension L/Q. 

Since < 7=1 (mod 4) we have that q splits completely in 
Q(*)/Q- Since q has order (p — 1) in Zp it follows that q 
remains inert in <Q{ojp)/<Q. Since q = 5 (mod 2®+^) and 5 has 
order 2 ® in Z 2 <i+ 2 , it follows that in the extension Q(w 2 '>+ 2 )/Q, 
q splits completely in Q(i)/Q but remains inert thereafter. 
Let q split in Q(*)/Q according to 

q = ttitt/ 

where tti = (a -f tb) and tt* = (a — ib). Now by using the 
fact that in a field tower [E : K : F] of field extensions, 

/e/f = Ze/k/k/f. 5 e/f = Se/kSk/f, [E : F] = /e/f pe/f, 

it follows that tti remains inert in the extension L/Q(z). 

To now find a non-norm element of unit magnitude, we note 
that since the units of Z[z] belong to the set {±1, ±*}, the 
associates of 

TTl = 


a lb belong to the set 















L 



Fig. 7. The inertial degrees for the maximal field of the CDA 


{a + lb, —a — ib, t(a + ib), —t{a + tb)}. 

It follows that since ab ^ 0, a — ib does not belong to the set 
of associates of a + ib. Our goal now is to show that 


TTl 



is a non-norm element, i.e., that the smallest exponent k for 
which 7 ^ is the norm of an element in L, is n. This is the 
case since if 

7 ^ = Ni^/r{i) some f G L 

then 

n—1 

1=0 

where tr is the generator of the cyclic Galois group of L/F. 
For i = f, a,b G Ol, we have, in terms of ideals of Ol, 

n—1 n—1 

1=0 1=0 

Since = tti we have that if (tti) divides {a^{x)) for 

some I and x G Ol, it must divide (a‘(x)), for all 1. This in 
turn implies that the power of (tti ) in the prime decomposition 
of (^i)"nr= _Q^(cr*(6)) is k mod n whereas the power of 
(tti) in the prime decomposition of =o^(CT'(a)) is a 

multiple of n. Equivalently k must be a multiple of n. 

When ni = 1, it is sufficient to take q to equal 5, and 
L = Q(a; 2 <i+ 2 ). The prime 5 splits in Q{i) as (1 -f 2i){l — 2i) 
and then each of (1 + 2 i) and (1 — 2 i) remain inert in the 
extension L/Q(^). The element 

1 + 2i 


is then a non-norm element for this extension, for the same 
reasons as above. This concludes the proof of Proposition Q 

□ 

b) Unit-magnitude, non-norm elements for F = Q(a; 3 ).‘ 
Let n = 2®ni, s G {0,1} where ni is odd. The proof is 
similar to when F = Q(*). Assume first that ni > 1. We 
hnd a prime p = 1 (mod ni), p > 3 and a prime q G Ij, 
<7=1 (mod 3), with ord(( 7 )|z ^ = ni. If s = 1, we also 
require that g = 3 (mod 4). Assume that we have found such 
a p and q. The argument for the rest of the statements in 
this paragraph are all exactly as in the case when F = Q(j): 
The conditions ord(g)|z * = ni and g = 3 (mod 4) (if s = 1) 


guarantee that the prime g remains inert in the ring of integers 
Ol' of the cyclotomic held L' = K(a; 2 s+i), where K is the 
unique subheld of degree ni in the extension Q(wp)/Q. The 
held U is cyclic over Q of degree 2'*ni. Since g = 1 (mod 3), 
the prime g splits into two distinct primes 7ri,7r2 in Z[a; 3 ] 
which are conjugates of each other. Let L = L'(w 3 ), which 
is cyclic over Q(a; 3 ) of degree 2®ni. Then tti will remain 
inert in the extension L/Q(w 3 ). The element 7 = ^ = -^ 
will then be a unit-magnitude (algebraic) non-norm element 
for the extension L/Q(a; 3 ), and the codes constructed with 
this data will then be full-diversity, full-rate, and have non¬ 
vanishing determinant, and of course, will satisfy the equal 
power-sharing constraint as 7 is of unit-magnitude. 

What is left is to hnd p and g. The prime p is found using 
Dirichlet as in the case where F = Q(t). To hnd g, hrst hnd 
an integer b that is simultaneously congruent to 1 (mod 3), to 
m (mod p), where m is a generator of , and (if s = 1) to 3 
(mod 4). This is possible by the Chinese Remainder Theorem. 
Next, hnd the prime g by applying Dirichlet’s theorem to the 
arithmetic sequence b-\-l{3p), I = 0 , 1 , 2 ,... if s = 0 and the 
sequence b -f ^12p), I = 0 , 1 , 2 ,... if s = 1. 

When ni = 1 (so s = 1), we take L to be Q(a; 3 )(z), and 
the prime g to be 7. Since g is congruent to 1 (mod 3) and 
to 3 (mod 4), g splits into 3 -f a ;3 and 3 -f in (^{ 103 ) but 
remains inert in the extension Q(i)/Q. It follows that each of 
3 -f u ;3 and 3 -f remain inert in the extension L/Q(u; 3 ). 
The element 

3 -I- a;3 

—2 

3 “h 

will then be a non-norm element for this extension, for the 
same reasons as above. This concludes the proof of Proposition 

E □ 


Appendix II 

Orthogonal Lattices in C>k, Where K/Q is Cyclic 
Galois oe Odd Degree 

We here show that the construction in [10] (of which a 
detailed exposition has been provided in [12, Section 5]) of 
lattices that belong in a cyclic Galois extension K of prime 
degree g over Q, actually gives without any modification 
orthogonal lattices for any odd degree n. We will follow the 
exposition in [12] closely, retaining even the notation in [12], 
and show that the proofs there only use the fact that n is odd, 
and not that it is an odd prime. 

To this end, let n > 3 be a given odd integer, and hx a 
prime p = 1 (mod n). Note that the existence of such a p 
is guaranteed since the sequence (1 + dn, d = 1, 2, • • • }, as 
shown by Dirichlet, contains infinitely many primes. Let w be 
a primitive p-th root of unity. Thus, Q(w) is cyclic of degree 
p —1 over Q, and contains the real sub field Q(a;-l-a;“^) which 
is cyclic of degree {p — l)/2 over Q. Since n divides p—1, 
there is a unique held K contained in Q(w) which is cyclic of 
degree n over Q. This is the held we will work with. Note that 
since n is odd, n divides (p— l)/2 as well, so K is contained 
in the real subheld Q(a; -f 0 ;“^). 

Recall that we are following the notation in [12]. Let G = 
Gal{Q{uj)/Q), with generator cr, chosen so that a{Lu) = w’’. 




where in turn, r is a generator of Z*. We let m = and 
observe that r™ = — 1 (mod p). We also choose A so that 
A(r — 1) = 1 (mod p). 

We dehne a by a = ~ The following result 

is just a combination of Lemmas 3 and 4 of [12], and since 
they have to do purely with the cyclotomic extension Q(a;)/Q 
and have nothing to do with n, their proofs remain valid; 
Lemma 8: The following equalities hold; 

1) (T(a) = 

2) a{uj^a) — uj^a 

3) = (-!)> 

We now dehne z = — to) G C>Q(tj), and 

(p-l)/n 

a: = TrQ(^)/K(z) = 

Note that x is in Ok, as z is in C2 q(cj)- Observing that 

GNGjj{iJ) = TrK/Qia\x)a\x)), 

we are interested in TvK/qixa*(x)). The following, which is 
Proposition 2 of [12], gives us the key to constructing the 
orthogonal lattice. 

Proposition 9: TrK/qixa*(x)) = p^So^t, for t = 0,..., n— 

1. 

Remark 1: Note that TrK/iQ{cr'^{x)a^ {x)) = 

TrK/q^{xa^~'^{x)). Thus, if we embed Ok in R" via 
a ^ v{a) = [a, (T(a),..., cr”“^(a)] (note that K is 

a real held), this Proposition says that the vectors 
[u(a;), u(cr(a:),..., u(cr"“^(a;))] are orthogonal to one 
another. 

Proof: For n being odd, we have 

n —1 

rrK/Q(a;cr*(a;)) = '^a°-{xa\x)) 

a—0 

n —1 (p—l)/n 
a=0 c,j = l 

and from Lemma IS] 

n-l {p-l)/n 

TrK/Q{xcr*{x)) = Y^ ^ 

a=0 c ,_7 = l 


Now the term (1 ~ 0 ;’’”^°") can be 

rewritten as (1 — 0 ;’’”*^°") l)a. Since n is 

odd, {p — l)/n is even, and hence, there are as many 

positive as negative terms in the expression 
and thus, the sum becomes zero. Similarly, the term 

becomes zero; this is because the terms in 

I]r=o(“ independent of 

c, while the term = 0 as {p — l)/n is even 

and there as many positive as negative terms. We thus hnd 

(p-l)/n 

TrK/^{xa\x)) = {-ir^p ^ (- 1 )=- 

C^l 

n-l (p-l)/n 
a=0 j=l 


We now have the following; 
Lemma 10: 


(P—!)/'”■ n — 1 (p—l)/n 

E E <-i) 


a+cn.a+t+jn 




c—1 a=0 j—1 

{p-l)/n n-l{p-l)/r~ 


= E (-d'E E 

fZ—1 a—0 k—1 

{p-l)/n n-l{p-l)/n 

= E <-i)‘'E E 

d—1 a—0 k—1 

Proof: See Appendix 11111 
As in [12], we write 


+ n<i+Tifc_j_^a + f + Tife 


UJ 




(p-l)/n n-l {p-l)/n 

E (-D'E E ‘ 






a=0 A;=l 


(p-l)/n (p-1) 

= E (-i*' E 


UJ 


d,t 


d^l 




where ujci,t = u;G"''^+P)^ and of course. 


We observe that since n is odd, (—1)°" = (—1)'^ and 
(-l)f" = (-l)J. Moreover, (-1)“(-1)“ = 1, and (-1)‘ 
is common to the sums above. By Lemma 0 we may replace 
by (—l)"*p. Thus we hnd, after rearranging the sums, 

that 


(p-l)/n 

TrK,^{xa\x)) = {-l)\-irp ^ (-!)=■ 

rn-1 (p-l)/n 
• ^ ^ (-1)^(1 
-a=0 j=l 
n-l (p-l)/rt 


-EE 

a=0 j=l 


^a + cn_^^a+t+jn 


(p-i) 


E<. 


p-1 if iUd^t = 1, 
— 1 otherwise 


To determine when t = 1, note that this happens (as in 
[12]) when t = nd—m + ki{p— 1). Since n is odd, n divides 
m, so n must divide t. This forces f = 0. 

We now have 0 ;^ * = 1 implies = — 1 (mod p), and 
writing —1 as r™, yields nd — m = l{p — 1) for some 1. This 
then gives d = [p — \){2l + l)/2n, which we may write as 
{21 -f 1) times {p— l)/2n) (note again that since n is odd, n 
divides (p—1)/2). Since d varies in the range 1,..., (p— l)/n, 
we hnd that I must be zero, that is, d = {p — l)/2n. Thus, 
= 1 precisely when t = 0 and d = {p — l)/2n. 




In particular, when t ^ 0 then 1 and we have that 

(p-l)/n (p-1) 

Tr^/Q{xa\x)) = (-l)‘+> ^ (-1)'^ ^ 

d—1 s—1 

{p-l)/n 

= (-!)*+> (-1)"(-1) 
d^l 

Once again, since n is odd, (p — l)/n is even, so the term 
^^^-i)/n(_^)d ^ t^O, TrK/Q{xa*{x)) = 0. 

When t = 0, we find 


(p-l)/n 


Tr]K/Q(a:cr*(x)) = (-l)’"p ^ 

d—l,d^{p—l) 1‘ln ^ 


(-i)'^(-i) 


and the right side then yields p + p(p — 1) = p^. To see this 
last fact, consider first the case where (p — l)/2 is even (i.e., 
p = 1 mod 4). Then, since n is odd, (p — l)/2n is also even. 
The sum equals - 

(_l)(p-i)/2™, Yve have already seen that, again because n 
is odd, Thus the right hand side in the 

equation above for TrK/Q(xcr*(x)) indeed yields p^ in this 
case. We can similarly deal with the case when (p — l)/2 
is odd (i.e., p = 3 mod 4), to find that in both cases, indeed 
T'fK/Q(a;cr‘(a::)) = p^ when f = 0. This proves the Proposition. 


Appendix III 
Proof of LEMMAfTol 

We wish to prove; 


P —1 1 P —1 

n —1 

c—1 a—0 j —1 


. + cn_^^a + t+jn 


P - 1 P — 1 

—;r~ n — 1 


d—1 a=0 fc=l 


Set m = and denote 'Ljm’L by In the above 

equation, the dependence on c, j, d, k is only through their 
values (mod m) or through their values (mod 2). If we assume 
2|m, which follows from the assumption that n is odd, we can 
then treat c,j,d,k as elements of Z^. We thus have 


c—1 a—0 j — 1 

n—1 

= E(-i) 

cGZm k^Zm 


a+cn . a+t+kn 


+r 


= E E<-')"‘E 


LJ 


j,a + cn_^^a+t+kn 


cGZm kGZ„ 


a—0 


We now make the change of variables: c = d + k (mod m) 
which implies, since 2|m, that c = d + k (mod 2) and hence 


d = c—k = c+ k (mod 2). As the pair (c, k) varies over all 
of (Zm X Zm), so does the pair {d, k). We thus have 


n — 1 


c€Zm kGZrr 


E E(-ir'‘E“ 

zGZm 

n—1 

= Ei-p'E E 


r +r 




d^Zm a=0 k^Zm. 

n—1 

deZm a=0 feGZm 


n—1 ■ 


“Et-D^EE' 

d—1 a=0 k—1 
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Appendix IV 

Proof of Approximate Universality and 
Information Losslessness of Perfect Codes 
(TheoremIUi 

The approximate universality part of the proof, is based on 
the derivation of the approximate universality conditions in 
[18]. It is reproduced here for completeness. 

Proof: Let Ai < A 2 < • • • < and > (2 > ■ ■ • > be 
the ordered eigenvalues of and respectively. 

Irrespective of the statistics of the channel, in the high-SNR 
regime, the probability of no-outage at multiplexing gain r, is 
shown in [3] to satisfy 

Pr(no-outage) = Pr ^^(1 + SKRA^) > ln(SNR’’)|, 

where n' := uim{n,nr). Through the Lagrange multiplier 
technique we determine 

n 

worst ~ / 2 

outage 

i—1 

by writing the functional as 

n n 

J(Ai,..., A„') = + P E/ lu(l + SNRAi) — /rrlnSNR 

i=l i=l 

and differentiating w.r.t. Ai, we obtain Ai = {n/h — SNR“^). 
We then use the Kuhn-Tucker conditions to verify that the so¬ 
lution \i = {^/li — SNR~^)+ is what gives the worst possible 
^1,worst’ for P such that 

n 

^ln(l-f SNR(p/Zi - SNR"^)+) =rlnSNR. 

i=l 

Solving the above, we obtain that 
G 

fx = SNR"(^"^) ^ If and thus Ai = ^ - 

i—1 

Substituting this value of Xi in d% and setting (i|; > 

SNR*^ for some e > 0, we obtain a condition on the smallest 
n' eigenvalues of the code m=ih > SNR” “E a condition 
satisfied by CDA codes with non-vanishing determinant. [16]. 






Now moving to perfect codes, we follow the approach in 
[22] to show that the linear dispersion matrices are unitary. 
This property, together with the full rate condition, establish 
the information losslessness and the entire theorem. 

As the code maps information elements, we consider 
n linear dispersion matrices of dimension n? x n. 

Starting with 

An,i = diag (G(*))„x„ , i = 1, • • • ,« 

where G{i) represents the row of G in we recursively 
create 

A,,, = u = l,2,---,n. 

Finally A„ is constructed as 

Au,o 


Au,n—1 

It is easy to see that the unitary nature of F and G makes each 
of the Au unitary 

A^^Au = In- 
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